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ABSTRACT 

The intergalactic neutral hydrogen which is responsible for the Lyct forest of quasar 
absorption lines is a tracer of much larger amounts of ionised hydrogen. The ionised 
component has yet to be detected directly, but is expected to scatter CMB photons via 
the Sunyaev-Zel’dovich effect. We use hydrodynamic simulations of a LambdaCDM 
universe to create mock quasar spectra and CMB sky maps from the same volume of 
space. We find that the high redshift Lya forest gas causes temperature fluctuations of 
the order of 1 /iK rms in the CMB on arcmin scales. The kinetic and thermal Sunyaev- 
Zel’dovich effects have a similar magnitude at redshift three, with the thermal effect 
becoming relatively weaker as expected at higher redshift. The CMB signal associated 
with lines of sight having HI column densities > is only marginally stronger 

than that for lower column density sightlines. There is a much more significant depen¬ 
dence of rms temperature fluctuation on mean Lya absorbed flux, however, suggesting 
that the CMB signal effectively arises in lower density material. We investigate the 
extent to which it is possible to cross-correlate information from the Lya forest and 
the microwave background to detect the Sunyaev-Zel’dovich effect at redshifts 2 — 4. 
In so doing, we are enable to set direct limits on the density of diffuse ionised inter¬ 
galactic baryons. We carry out a preliminary comparison at a mean redshift 2 = 3 
of 3488 quasar spectra from SDSS Data Release 3 and the WMAP first year data. 
Assuming that the baryons are clustered as in a LambdaCDM cosmology, and have 
the same mean temperature, the cross-correlation yields a weak limit on the cosmic 
density of ionised baryons Dbj. As a fraction of the critical density, we find Hb.i < 0.8 
at 95% confidence. With data from upcoming CMB telescopes, we anticipate that a 
direct detection of the high redshift ionised IGM will soon be possible, providing an 
important consistency check on cosmological models. 
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1 INTRODUCTION 

The Lya forest seen in quasar spectra is now thought to 
be produced by neutral hydrogen atoms in an intergalac¬ 
tic medium (IGM) which is in photoionisation equilibrium 
(see e.g., Rauch 1998 for a review). At low to moderate red- 
shifts, the fraction of neutral atoms is small (1 in ~ 10®; e.g. 
Gunn & Peterson 1965) but knowledge of the photoionisa¬ 
tion rate (from summing the contribution of known sources; 
e.g. Haardt & Madau 1996) makes it possible to extrap¬ 
olate from this one part in a million measurement to an 
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estimate of the total baryon density Q.b (Weinberg et al. 
1997, Rauch et al. 1997). Loeb (1996) has argued that the 
ionised component of the IGM should be directly observ¬ 
able by inverse Gompton or Doppler scattering of cosmic 
microwave background (GMB) photons by free electrons 
(Sunyaev-Zel’dovich [1972,1980]; hereafter, the SZ effect). 
Such a detection would provide consistency checks on the 
currently favoured cosmological model for structure forma¬ 
tion, and place constraints on the IGM density, tempera¬ 
ture, ionisation state and velocity field. In this paper, we 
use a hydrodynamic cosmological simulation to investigate 
how observations of the CMB temperature and Lya spectra 
of quasars drawn from the same region of sky can be used in 


2 R.A.C. Croft ei al. 


the search for the hitherto unseen 99.9999% of the hydrogen 
associated with the Lya forest gas at redshifts 2 « 2 — 5. 

In favoured models for structure formation, the optical 
depth for Lya forest absorption arises in a continuously fluc¬ 
tuating medium (e.g., Cen et al. 1994; Zhang, Anninos, & 
Norman 1995; Petitjean, Miicket, & Kates 1995; Hernquist 
et al. 1996; Katz et al. 1996a; Wadsley & Bond 1997; Theuns 
et al. 1998; Dave et al. 1999). The density of neutral hydro¬ 
gen is inversely proportional to the photoionisation rate P: 


r = 
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and the optical depth to absorption is given by 

TLya = - f \H {z)n, (2) 

rrieC 


where / = 0.416 is the Lya oscillator strength, n is the neu¬ 
tral hydrogen number density, and A = 1216 A (Gunn and 
Peterson 1965). In order for r be of order unity, (as seen, for 
example at redshifts 2 ~ 3), n must be much smaller than 
the total density of hydrogen atoms, as noted earlier. The 
gas in photoionisation equilibrium lies at densities ~ 1 — 10 
times the mean, and obeys a power law relationship between 
density and temperature (Gnedin & Hui 1998). At 2 = 3, 
according to hydrodynamic simulations (Dave et al. 2001), 
more than 90% of the baryons are expected to be in this 
diffuse component. The tiny neutral fraction is readily de¬ 
tectable, but not the dominant ionised portion of the mate¬ 
rial. Owing to gravitational evolution, at 2 < 2 most of the 
baryons are expected to reside in a shock heated IGM, with 
little neutral hydrogen (Cen & Ostriker 1999, Dave et al. 
2001). 

The interaction of CMB photons with free electrons in 
the intergalactic plasma was first considered by Sunyaev and 
Zel’dovich (1972). Inverse Compton scattering preferentially 
increases the energy of CMB photons, while conserving pho¬ 
ton number, leading to a spectral distortion whose ampli¬ 
tude is proportional to the product of electron temperature 
and density (the thermal SZ effect). Doppler scattering in¬ 
duces an intensity fluctuation with the same spectral shape 
as the CMB itself (the kinetic SZ effect). The thermal effect 
is one of the main sources of secondary CMB anisotropies 
on small angular scales. 

The perturbation in the CMB thermodynamic temper¬ 
ature resulting from scattering of nonrelativistic electrons 
is 


AT 
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~ —2y for X « 1, 
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where x = hu/kTcms — iz/56.85 GHz is the dimensionless 
frequency and the second expression is valid in the Rayleigh- 
Jeans limit, which we assume henceforth. The quantity y is 
known as the comptonization parameter and is given by 

f ,.nek{Te — Tcmb) 

y= / dl - 2 -> (3) 

where the integral is performed along the photon path. 

The kinetic SZ effect arises from the motion of ionised 
gas with respect to the rest frame of the CMB. The resulting 
temperature fluctuation is ATjT — —b where 



gives the magnitude of the effect and Vr is the component of 
the gas peculiar velocity along the line of sight (positive sign 
for receding gas, negative for approaching) to the observer. 

The gas which causes Lya forest absorption has a den¬ 
sity near the cosmic mean, with a volume weighted temper¬ 
ature around 2 x 10'* K at 2 = 3 (e.g, Schaye et al. 2000, 
McDonald et al. 2001). At this temperature, we expect the 
contribution from the kSZ to dominate the rms tempera¬ 
ture fluctuations of the CMB from the Lya forest. Using 
observed Lya forest line counts, Loeb (1996) estimated that 
the rms temperature fluctuations from the kSZ would be 
AT/T ~ 10“® for the Lya forest lying between 2 = 2 and 
2 = 5. This is for angular scales of order 1 arcmin, with 
a decrease on larger scales. There is much competition 
from other signals on these small scales (including low red- 
shift thermal SZ from galaxy clusters, the Ostriker-Vishniac 
(1986) effect, patchy reionisation (see, e.g. McQuinn et al. 
2005), and foreground sources). The suggestion to cross¬ 
correlate information from SDSS Lya forest spectra and the 
CMB was also made by Loeb (1996) in order to try to better 
extract the signal. In this paper, we measure the kSZ effect in 
simulations, as well as the thermal effect. The mass weighted 
temperature of the IGM at 2 ~ 3 is expected to be closer to 
IQ^K (e.g., Springel, White and Hernquist 2001), and so the 
thermal effect should be significant. How much of that hot 
gas is physically associated with Lya absorption will affect 
how well the tSZ can be detected by cross-correlation. 

We note that both being absorption phenomena, neither 
the SZ signal nor the Lya forest opacity are directly affected 
by the inverse square law. The problem of finding bright 
QSO background sources aside, they therefore both have an 
advantage in the hunt for baryons at high redshifts. Indeed, 
the value of the SZ effect for finding galaxy clusters at the 
highest redshift has long been recognised (see e.g., Carlstrom 
et al. 2002 for a review). We aim to explore its potential for 
finding diffuse ionised material at 2 > 2. 

The plan of this paper is as follows. In Section 2, we 
give details of the hydrodynamic cosmological simulation 
and how mock Lya spectra and CMB temperature sky maps 
were constructed from the same volume of space. In Section 
3, we compute statistical measures of the absorption and 
examine correlations between the two probes of the inter¬ 
galactic medium. In Section 4, we compare directly to ob¬ 
servational data from the SDSS and WMAP, and in Section 
5 we summarise and discuss our results. 


2 SIMULATED SPECTRA AND SKYMAPS 
2.1 Hydrodynamic simulation 

We employ a hydrodynamical simulation of the popular 
ACDM cosmology to make our prediction of the Lya and 
the SZ effects. In order to both resolve the Lya features 
and the bulk velocities responsible for the kSZ effect, we 
require a large simulation volume with good mass resolu¬ 
tion. The data snapshots we use are from the simulation 
“G6”, also used by Nagamine et al. (2005) to which we 
refer the reader for more details, which describes a repre¬ 
sentative cubic volume of space 100 h,“*Mpc on a side, in 
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a cosmological constant-dominated CDM nniverse, consis¬ 
tent with measurements from the WMAP satellite (Bennett 
et al. 2003, Spergel et al. 2003). The relevant parameters 
are Qa = 0.7, flm = 0.3 fib = 0.04, and a Hubble constant 
Ho — 70 km s“^ Mpc“^. The initial linear power spectrum 
is cluster-normalised with a linearly extrapolated amplitude 
of as = 0.9 at 2 = 0. The simulation was performed with 
the smoothed particle hydrodynamics (SPH) code Gadget- 
2 (Springel et al. 2001, Springel 2005), which is based on 
the entropy-conserving approach of Springel & Hernquist 
(2002), and is an extension of the “G-series” runs of Springel 
& Hernquist (2003b). 

The baryonic matter is initially represented by 486® 
SPH particles and the collisionless dark matter by 486® N- 
body particles. The initial mass per gas particle is therefore 
9.7 X IO^Mq and dark matter 6.3 x IO^Mq. The gravita¬ 
tional softening length is 5h~^kpc. The calculation includes 
gas dynamics, cooling, a multiphase treatment for star for¬ 
mation (see Springel & Hernquist 2003a), and a uniform 
background of ionising radiation in the manner of Katz et 
al. (1996b), renormalising the Haardt & Madau (1996) spec¬ 
trum to reproduce the observed Lya forest mean optical 
depth. A study of the SZ effect using this code was carried 
out by Springel, White, & Hernquist (2002), who noted that 
the additional physics beyond adiabatic gas dynamics had 
a relatively small effect; changing, for example, the mean 
comptonization by 20%. 

2.2 Mock Lya spectra 

We use simulation outputs from redshifts 2 = 2, 3, 4, 5 and 6 
to make mock Lya spectra. For each output, this was done 
in the usual manner, by integrating through the SPH kernels 
of the particles to obtain the neutral hydrogen density held, 
and then convolving with the line of sight velocity held (see 
e.g., Hernquist et al. 1996). At each redshift, we compute 64® 
lines of sight, equally spaced on a square grid. In order to ap¬ 
proximate longer spectra, we join together spectra after per¬ 
forming a random translation, rehection and choice of axis, 
as is done for example when making skymaps by ray-tracing 
through multiple copies of the simulation volume (e.g., da 
Silva et al. 2000). As our hducial spectrum length, we take 
the distance between the Lya and Ly/3 emission lines. Later 
in the paper, when we cross-correlate the Lya forest and 
CMB at a particular redshift, our analysis will refer to this 
Lya to Ly/3 region, centered on the redshift in question. For 
reasons of computational simplicity, we choose to ignore evo¬ 
lution over the length of each spectrum, both in clustering 
and in the mean optical depth. The Lya to Ly/3 regions cor¬ 
respond to comoving spatial distances of 343, 343, 330 and 
314 /i“^Mpc when centered on redshifts 2 = 2,3,4 and 5, 
respectively. We add partial spectra to the full spectra to 
ensure the correct overall length. 

For each spectrum, we compute the mean effective op¬ 
tical depth over its entire length, i.e., 

reff = -ln(F’), (7) 

where (F) is the mean absorbed flux in the spectrum, 
Throughout most of this paper, we will use this quantity, 
Teff, to determine the relative amount of Lya absorption in 
each sightline. Another possible measure is to compute the 
total column density of neutral hydrogen along each sightline 


from Lya to Ly/3 . Being unaffected by saturation, this will 
tend to reflect the presence of individual high density clouds. 

As noted earlier, the simulation was run using the UV 
ionising background of Haardt & Madau (1996). In our 
subsequent analysis, we rescale the Lya optical depths by 
a small amount after computing the spectra by adjust¬ 
ing them slightly so that {F) averaged over all spectra re¬ 
produce the observational results of Schaye et al. (2004), 
{F) = 0.878, 0.696, 0.447 for 2 = 2, 3, 4. 

2.3 CMB skymaps 

At each redshift, we construct maps of the Gompton y pa¬ 
rameter and the Doppler b parameter by projecting the ap¬ 
propriate SPH kernel weighted quantities along one axis. 
Unlike Springel et al. (2001), for example, we do not project 
along lines which meet at a single 2 = 0 observer, as we 
will be interested only in the SZ signal generated by gas 
which falls within the same redshift interval as our Lya for¬ 
est spectra. In order to compare GMB maps and spectra 
for the same volume of space, we randomly translate and 
reflect our SZ maps in the same manner so they match up 
with the spectra from the previous section. We generate 40 
maps, which have an angular size which subtends the whole 
box at each particular redshift; i.e., 94.7 arcmin at redshift 
2 = 2, and 77.2, 68.5, 63.1 and 59.5 arcmin at 2 = 3, 4, 5 and 
6 respectively. We make maps at a resolution corresponding 
to 2048® pixels, but also average the pixels so that we have 
64®, the same number as the Lya spectra. 

2.4 Examples 

Looking at kSZ maps and Lya spectra generated from the 
same region of space, we can see if there are any obvious 
visual relationships between the GMB temperature in pixels 
and the appearance of the absorption features. In Figure 0 
we show the Gompton y parameter, from the thermal SZ 
effect at 2 = 2, and alongside it five spectra which are taken 
with the sightline along the axis perpendicular to the plane 
of the map. Here, in both maps and spectra, we do not plot 
the entire Lya to Ly/3 region (we do this later, below), but 
just the material in one single cubic simulation volume. It is 
apparent from the map that much of the volume contributes 
to a very small y parameter, and that most of this redshift 
slice has y < 10“® or below. There are also a number of 
obvious galaxy clusters, with y > 10“® 

Of course, a map of the y parameter from all material 
between 2 = 0 and the last scattering surface would have 
quite a different appearance, with the mean y parameter 
being higher, of the order of 10“®, and with many of the 
obvious filaments seen in Figure^washed out by projection 
effects. Figure 2 of White et al. (2002) demonstrates this, 
using a similar simulation. 

Sightlines (a)-(c) in Figure 0 were chosen from those 
among the 64® we made which have large values of y. Spec¬ 
trum (a) is visibly taken from near a cluster, seen in the SZ 
map. There, however, does not seem to be much evidence 
of this in Lya absorption, at least in terms of saturated ab¬ 
sorption. The other two spectra, (b) and (c) are also fairly 
similar, even though they are both associated with large y 
values. Because there is not much saturated absorption, we 
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Figure 1. The tSZ CMB y parameter (left panel) for a projection of one simulation volume (side length 100 h“^Mpc). We show the 
positions of five pixels for which we have drawn Lyo sightlines (into the page). The five spectra are shown in panels (a-e). Black curves 
in the right panels are the mock transmission spectra, while the dashed red curves are the peculiar velocity. Sightlines (a) and (b) and 
(c) have been chosen so that y is large (for (a), y = 5 X 10~®, for (b), y = 1.4 X 10“®), and for (c), y = 3.3 X 10“®), all among the top 
1 percent of y-values. Sightlines (d) and (e) have small values ol y, y = 5.13 X 10“^® and 3.8 X 10“^® respectively. Sightline (b) is the 
same as sightline (b) in Figure|^The Lya spectra are of length 9900 km s“^ (one simulation box length). 
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Figure 2. Projection through the simulation box showing the kSZ CMB Doppler parameter, b (left panel). Pixels with a positive values 
of b contain material moving on average away from the observer, and negative b values towards the observer. We show the positions of 
five pixels for which we have drawn Lya sightlines (into the page). The five spectra are shown in panels (a-e). Black curves in the right 
panels are the mock transmission spectra, while the dashed red curves are the peculiar velocity. Sightlines (a) and (b) have been chosen 
so that b is large and positive (for (a), b = 7 X 10“® and for (b), 5 = 4 X 10“®), both among the top 1 percent of b values. Sightline 
(c) is the pixel with the smallest value of |6|, —8 X 10“^^. Sightlines (d) and (e) have large -ve values of b, —3 X 10“® and —6 X 10“® 
respectively. The Lya spectra are of length 9900 km s“^ (one simulation box length). 



can see that it is unlikely that column density will correlate 
well with SZ decrement. 

Sightlines (d) and (e) in Figure 0 have been chosen to 
have values of y near the bottom of the rank-ordered list of 
all pixels. From the SZ plot it can be seen that the spectra 
lie on sightlines though obvious voids, and the y values are 
10000 times smaller than for panels (a)-(c). On examina¬ 
tion, these spectra appear to have very few wide absorption 
features compared to spectra (a)-(c), and the overall level of 
Lya absorption is slightly lower (particularly (e)). 

Turning now to the map of the kSZ effect (Figure |2|l, 
we see that much of this redshift z = 2 slice has a Doppler 
|6| value of ~ 10“®, although there are several clusters and 
filaments dense enough and moving fast enough to reach 
l&l = 10“®. Clusters can be seen moving towards the ob¬ 
server (blue) and away (red) . The large coherence length of 


the velocity field can be seen also, from the size of the red 
and blue structures. Projecting the kSZ effect for a larger 
volume, for example subtending z = 2 — 5, will lead to more 
structures but also more cancellation of positive and nega¬ 
tive b values. The general level of fluctuations likely agrees 
at least roughly with the estimate of Loeb (1996) {rms 
6 ~ 10 “®). 

Again, the kSZ plot can be compared to full lightcone 
maps (e.g. Figure 2 of Springel et al. 2001), which show a 
similar pattern of approaching and receding material, also 
with less obvious filamentary structure. Sightlines (a) and 
(b) in Figure 121 were chosen to have large positive Doppler 
b values, and (d) and (e) negative. It can be seen that their 
spectra contain fairly substantial amounts of absorption, 
and one can understand how the mass weighted velocity will 
give rise to the b values seen. Panel (c) has a velocity field 
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which cancels in the left and right halves of the spectrum, 
yielding a b parameter much closer to zero (b ~ —10“^^). 

The number of sightlines passing close to galaxy groups 
and clusters in Figures 0 andis relatively small, at least in 
this redshift interval. This means that any clustering statis¬ 
tic which weights spectra equally will tend to be sensitive 
to physical conditions in the regions with relatively small 
overdensity. 

Looking at a few individual sightlines, it seems likely 
that there is a connection between Lya forest absorption 
and SZ-induced temperature fluctuations. One way to ex¬ 
amine this further would be to look at skymaps of both the 
Lya forest and CMB. Of course, with real observations it un¬ 
likely that we would have a dense enough grid of background 
quasars to make a true Lya map, but with the simulations, 
we can use such a technique to examine the spatial pattern 
of absorption. In Figure |21 we show examples of such maps, 
plotting the thermal and kinetic SZ effects alongside angular 
maps of the Lya optical depth and column density. It should 
be noted that we plot the absolute value of the temperature 
change produced by the kinetic effect. Unlike the previous 
two figures, here we show a projection through the volume 
of space corresponding to an entire Lya to Ly/3 region (not 
just one simulation box). The results are centered on red- 
shift z = 2, and the angular size of the maps is 95 arcmins. 
Because we have a grid of 64^ Lya spectra, we have rebinned 
the CMB maps so that they have the same resolution (1.5 
arcmin pixels). 

The grayscales of the SZ maps are the same, and so we 
can see from the left two panels of Figure|^that the thermal 
and kinetic effects trace similar features, with the thermal 
effect being more prominent at this redshift. Because of the 
differing frequency dependence of the two SZ effects, cross¬ 
correlation of the two maps may eventually be possible us¬ 
ing observational data, as suggested for example by da Silva 
et al. (2000). In the present context, we are interested in 
the potential cross-correlation of the Lya forest spectra and 
the SZ effects. In the third panel of Figure El we show the 
mean effective optical depth, given by Equation o for spec¬ 
tra through the center of each pixel. It can clearly be seen 
that the patches of absorption largely follow the SZ tem¬ 
perature well, although there are several regions with large 
differences. The Lya Tefr fluctuations have a somewhat lower 
contrast, varying by about a factor of two across the map. 
The obvious groups and clusters in the SZ maps do not show 
up as prominently in the Lya map. This is likely to at least 
partly owe to the fact that regions with large amounts of 
Lya absorption are saturated. 

This is not the whole story, however, as we can see by 
looking at the rightmost panel of Figure El Here, we show 
the HI column density in each spectrum integrated from 
Lya to Ly/3 . There is much smaller scale structure present, 
as the dense knots of neutral hydrogen are picked out. The 
range of column densities varies from Nhi = 10^"^ cm“^ 
to 10^® cm“^. The structures in Nhi do not correlate well 
with the CMB maps. The smooth, large-scale structures in 
the SZ effect maps (and Lya Tea ) are not readily seen. 
This observation will guide our use of Tefi as the main Lya 
statistic to be used. 

In Figure El we show a scatter plot of the Lya Tea 
against the thermal SZ temperature decrement, at redshifts 
2,3 and 4, again from maps with 64^ pixels (1.5, 1.2 and 



Figure 4. Thermal SZ temperature decrement Ttsz vs Lya 
forest Tea for spectra at three different redshifts. 

1.1 arcmin size for z=2,3,4). The correlation between the 
two quantities is apparent, and in particular at each redshift 
there are few lines of sight with very low values of Tefi but 
substantial SZ decrements. There is a large amount of scat¬ 
ter, however, indicating that even if foreground/background 
noise is not an issue, it will be necessary to average over 
many sightlines to obtain statistically significant results. 


3 STATISTICAL MEASURES 

In this section, we use the set of skymaps and spectra gener¬ 
ated from our simulation to examine statistically the CMB 
fluctuations caused by Lya forest gas. 

3.1 Rms temperature fluctuations 

In Figure El we show the rms temperature fluctuations av¬ 
eraged in top-hat circles on the plane of the sky, calculated 
from maps such as those shown in the left two panels of Fig¬ 
ure El We show results centered on integer redshifts between 
z = 2 and 2 : = 6, in each case, calculating only the CMB 
fluctuations caused by gas within the Lya to Ly/3 region of 
the spectrum. The thermal and kinetic SZ effects are shown 
separately, and it can be seen that at 2 = 2, Trma is approx¬ 
imately three times larger for the tSZ effect, as we expect 
from Figure El At 9 = 1 arcmin, the rms T fluctuations are 
~ 2^K for the tSZ at this redshift. 

The curves for all redshifts have an angular dependence 
given by ~ At 2 = 3, the two SZ effects are compa¬ 
rable, and beyond that the kSZ effect stays approximately 
constant out to 2 = 5. The rms peculiar velocity of the 
gas increases slightly with decreasing redshift, as detailed 
in Tabled This is balanced by the decrease in angular size 
for a given comoving length which shifts the curves to the 
left. The main change in the kSZ at 2 = 6 is caused by the 
switching on of the UV ionising background in our model, 
and so the ionised fraction is small then. The temperature of 
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Figure 3. Sky maps of CMB temperature (left panel) owing to the kSZ effect at z=2 (for a volume equivalent to the entire Ly« to Ly/3 
spectral width). The width of the panel is 95 arcmin. We also show |r|, mean flux < F > in Lya spectra and total column density in 
Ly« spectra. In each map, there are 64 X 64 pixels. 



9 (arcmin) 

Figure 5. Rms CMB temperature fluctuations in top-hat circles 
as a function of radius (x-axis) for three different redshifts. In each 
case, we plot the rms fluctuations associated with a voiume which 
wouid contain material with HI absorption waveiengths from Lyo 
to Ly/3 centered on the indicated redshift. Thin lines are kinetic 
SZ and thick iines thermal SZ. 

the gas influences the fluctuations owing to the thermal ef¬ 
fect directly, and this has a much more rapid dependence on 
redshift. In Tabled we can see that at z = 6, the electron- 
weighted mean temperature is > ICK, as the free electrons 
are concentrated in the few collisionally ionised structures 
present at that redshift.^ By 2 = 5 when most of the IGM is 
photoionised, the electron-weighted temperature is 2.6 x 
K, and as gravitational evolution continues, shock heating 

^ We ignore any difference between the electron and ion temper¬ 
atures; see Yoshida et al. (2005). 


Table 1. The mass density and eiectron number weighted gas 
temperature and velocity dispersion in the simulation at different 
redshifts. 


2 Tp Te '^rms,e 

(K) (K) (kms“l) (kms“l) 


6 

7.1 : 

X 

10® 

1.9 

X 

10® 

136 

166 

5 

2.57 

X 

10^ 

2.60 

X 

10^ 

149 

148 

4 

6.86 

X 

10"* 

6.89 

X 

10"^ 

164 

164 

3 

1.77 

X 

10® 

1.76 

X 

10® 

185 

185 

2 

5.04 

X 

10® 

5.04 

X 

10® 

214 

214 


increases the temperature by another factor of 20 by 2 = 2. 
The rise in T^ms owes mostly to this, with a small part com¬ 
ing from the factor 2.3 increase in linear density fluctuations 
themselves between the two redshifts. 

We note that although the thermal SZ effect is domi¬ 
nant at redshifts 2 < 3, the signal that is actually generated 
in gas that directly produces the Lya forest is a small frac¬ 
tion of the total thermal SZ signal at that redshift. This 
is because the temperature of Lya forest-producing gas at 
the mean density at 2 = 3 is ~ lO'^K, low enough that 
the thermal SZ effect should be extremely small (e.g., Loeb 
1996). The thermal SZ signal we are instead seeing is largely 
generated in the Warm-Hot Intergalactic Medium (WHIM, 
e.g. Dave et al. 2001) at these redshifts which is collisionally 
ionised, and is responsible for the higher mass-weighted tem¬ 
perature seen in Table 1. As we discuss below, the thermal 
SZ signal at this redshift is strongly correlated with the Lya 
absorption signal, indicating that the WHIM signatures are 
physically associated with the same large-scale structures as 
the Lya forest, tracing their higher density regions. 

3.2 Cross-correlation function 

We have seen that IGM gas at redshifts 2 = 2 to 2 = 5 
yields rms GMB temperature fluctuations on the order of a 
few pK. One could calculate the angular power spectrum, 
the Cis, and then see whether it could be separated from 
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Figure 6. Cross-correlation of |T| (for kSZ) or T (for tSZ) and 
F/ < F > as a function of pixel-pixel distance, for three different 
redshifts. Results for tSZ are shown by thick lines and kSZ by 
thin lines. 


other small angle fluctuations. Owing to the multitude of 
other possible signals, however, it is likely to be more useful 
to cross-correlate with material at a known redshift, and use 
this to extract information. As the Lya forest contains most 
of the baryons at these redshifts, and the mean effective 
Lya optical depth appears to correlate visually with CMB 
temperature in Figure |21 we use this quantity in our cross¬ 
correlation. 

We compute the cross-correlation separately for the tSZ 
and kSZ effects. For the kSZ, because the magnitude of the 
effect can be positive of negative, we work with \Ttsz\- Al¬ 
though this is not easily done with observational data, in 
the current section we are interested in making a relative 
comparison between the magnitude of the tSZ and kSZ ef¬ 
fects, where this is a useful approach. In order to do the 
cross-correlation, for each Lya spectrum we first calculate 
Areff = -(refi - (reff))/(reft), where (refi) = - In(F’) and we 
use the average mean transmitted flux for all spectra. The 
cross-correlation we compute is (Areff AT), which has units 
of pK. For the kSZ we use AT = |r| — (|T|) and for the tSZ, 
AT = T — (T) = T. The results for redshifts z = 2 — 4 are 
shown in Figure]^ We have not rebinned the CMB tempera¬ 
tures to the angular resolution of the spectrum grid, so that 
we are able to show results down to angular scales below 1 
arcmin. 

Unlike the case for the rms T fluctuations, one can see 
that the cross-correlation statistic for the kSZ effect depends 
on redshift, driven by evolution in the Lya forest. On first 
inspection, this is somewhat puzzling, as Lya Tefi fluctua¬ 
tions should be increasing towards high redshift, as the mean 
optical depth increases, but instead they are smaller. This 
is because we have divided out (reft), and so the fractional 
fluctuations are indeed smaller. Measurements based on the 
flux F, instead divide out {F) which is smaller at high red¬ 
shift (see e.g., Croft et al. 1998). 

This statistic and its spatial dependence could in prin¬ 
ciple be measured observationally for the thermal SZ effect. 
We will return to this point in Section 5. 


3.3 One dimensional distributions 

We have seen from Figure0]that there is a dependence of SZ 
decrement and the Lya Teff for spectra in the same pixel. We 
will average this relationship over many pixels and bin as a 
function of , in order to visualise the mean trend of T vs 
Teff . This simple statistic is one which we will concentrate on 
in this subsection, and use in Section 4 below to compare to 
observational data. The mean trend of T vs Teff does depend 
on the resolution of the CMB map, and we present results 
for the same resolution as the spectrum grid (64^ pixels, 
each of width 1.2 arcmin at z=3) as well as smoothing with 
a Gaussian filter with FWHM 13 arcmin (the WMAP W- 
band resolution, see Bennett et al. 2003). We will also plot 
the rms temperature fluctuations in pixels caused by the 
kinetic SZ effect. 

Figure | 7 | shows our results for three different redshifts, 
with the top panels being the tSZ decrement. We can 
see that the overall level of the temperature decrement is 
strongest for z = 2 and weakest for 2 = 4, as expected. 
At z = 2, pixels with Lya Teff =0.07 have tSZ temperature 
decrements of ~ 0/rK, and the pixels with Teff =0.2 have 
decrements ~ —1.5/rK. Between these two extremes, the 
trend of T with Teff is close to linear. The Lya spectra with 
the very highest levels of absorption, Teff 0.2, exhibit a 
weakening of the decrement, likely owing to radiative cooling 
of the gas in these high density regions producing relatively 
more neutral hydrogen. The signal smoothed at the WMAP 
resolution is somewhat weaker, but still easily measurable 
from the simulations. 

In the middle panels of Figure |7| we show the kSZ in¬ 
duced rms temperature fluctuations. They exhibit a rise 
with increasing Teff , in accord with the rise in density probed 
by the sightlines, again with a turndown for the spectra with 
the very highest absorption levels. As with the rms in circles 
plotted in Figure the level of fluctuations changes little 
with redshift. Smoothing the map to the WMAP resolution 
does affect the level strongly, however, reducing the magni¬ 
tude of the signal by roughly a factor of ten. As rms fluctu¬ 
ations are already more difficult to measure in the presence 
of noise than the mean level, this means that it is likely that 
the thermal rather the kinetic SZ effect will be easiest to 
target when considering the analysis of observational data. 

One can ask whether the signal from the tSZ and kSZ 
effects arises in low or high density gas, and whether we are 
really sampling the low density IGM probed by the Lya for¬ 
est. Because our measure of the Teff - T correlation is angu¬ 
lar pixel weighted (we are using a histogram of pixel values), 
the fact that most of the sightlines we use will contain rel¬ 
atively little gas comes into play. For example, only a few 
of the pixels in Figure O pass close to clusters and groups, 
and most of the spectra with the lowest 90 percent of Teff 
values are likely to pass through only void-like regions in the 
Lya to LyP interval. This means that while the thermal SZ 
does directly probe the mass weighted (or rather electron 
number-weighted) IGM temperature, and so is dominated 
by virialised objects, by breaking up the signal as a function 
of Teff we introduce sensitivity to what is occurring in low 
density regions. 

In Figure 13 we saw that there appeared to be a much 
better correspondence between GMB temperature and Teff 
than HI column density. Nevertheless, from looking at the 
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Figure 7. Top panels: Dependence of mean temperature from the thermal SZ effect on Tgff = — log < F > for spectra running 
through the pixels. Middle panels: Dependence of rms temperature fluctuation in pixels (1.5, 1.2 and 1.1 arcmin size for 2=2,3,4) on 
Teff = — log < F > for spectra running through the pixels, for three different redshifts (left to right). The round symbols which make up 
the upper curve in each case show results for pixels, unsmoothed, and the triangles for maps smoothed with a Gaussian filter of FWHM 
13 arcmins, similar to the WMAP W-band resolution. Error bars are from the scatter between 10 fields in each case (size 95, 77, 69 
arcmin for z = 2, 3,4). Bottom panels: the pdf of Tgfj = — log < F > for spectra from Lya to Ly/3. 





Figure 8. Top panels: Dependence of rms temperature fluctuation in pixels (same size as Figure 7) on NHI for spectra running through 
the pixels, for three different redshifts (left to right). The round symbols which make up the upper curve in each case show results for 
pixels, unsmoothed, and the triangles for maps smoothed with a Gaussian filter of FWHM 13 arcmins, similar to the resolution of the 
WMAP satellite. Error bars are from the scatter between 10 fields in each case (size 95, 77, 69 arcmin for 2 = 2,3,4 ). Bottom panels: 
the pdf of NHI for spectra from Lyo to Ly/3. 


maps it was not immediately obvious that the very highest 
density objects are in fact associated with large amounts 
of HI and have a significant SZ signal. In Figure |H1 we plot 
CMB T vs Nhi, and see that there is a slight increase in 
the SZ signal for sightlines with Nhi = — 10*^® cm“^, 

for all redshifts. The curves, however, are relatively flat at 


higher column densities. It appears that information from 
saturated regions will not be helpful in carrying out Lya 
forest CMB correlations, and so there is little need for high 
resolution spectra. In order to shed further light on this, in 
Figurel^we have plotted values of Nhi vs Tea for individual 
sightlines, as a scatter plot (as always, we average over the 
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Figure 9. A^hi vs reff for spectra at three different redshifts. 

Lya to Ly/3 region in each spectrum). We expect there to 
be an envelope, a minimum Nhi for sightlines with a given 
Teff . This minimum level would correspond to a uniform 
distribution of absorbing material (see e.g., Weinberg et al. 
1997). We can see such an envelope, but also a large scatter 
towards higher values of Nhi- The fractional scatter in Nhi 
increases for larger Tefr , so it is not surprising that for high 
Teff there is little evidence of a trend in Figure |H] 

4 LIMITS FROM SDSS AND WMAP DATA 
4.1 Method 

The microKelvin fluctuations produced by the Lya forest 
gas are likely too small to be detected directly in current 
data (e.g., Bennett et al. 2003, Kuo et al. 2004). However, it 
is possible that by comparing observational statistics with 
simulation results we can come up with interesting limits 
on the population of ionised baryons, as well as estimating 
the properties of the observations necessary to make a di¬ 
rect detection. Based on our investigations in the previous 
sections, a good strategy will be to measure Teg for the Lya 
forest as well as the CMB temperature in pixels where there 
are Lya spectra. For the present limited observational study, 
we choose to investigate the one-dimensional distributions of 
these statistics, as in Section We do this for reasons of 
simplicity, and because breaking up the signal into sightlines 
with high and low absorption levels will enable us to look 
for an obvious distinct signature (unlike a full spatial cross¬ 
correlation analysis which would have more discriminating 
power, but be more difficult to interpret). 

In our analysis, we calculate the mean and rms CMB 
temperature for pixels binned using the Tes of the spectra 
that pass through them. We use a jackknife estimator (e.g., 
Bradley 1982) to compute the error bars, breaking up the 
data into 50 subsamples. We have also computed error bars 
by randomly rotating the CMB data with respect to the Lya 
forest data and find them to be broadly consistent. 
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4.2 Lya forest data 

We use the Lya forest of quasar spectra taken from Data 
Release 3 of the Sloan Digital Sky Survey (Abazajian et al. 
2005). Our data sample consists of 3488 quasar spectra with 
quasar redshifts greater than 2.4, pruned to remove Broad 
Absorption Line quasars. The particular dataset was com¬ 
piled from the publically available data by Scott Buries and 
kindly made available to us (Buries, private communica¬ 
tion) . Owing to the relatively high low-redshift limit for the 
quasars, the spectra consist of the full Lya to Ly/3 region of 
the forest, lying above the atmospheric cutoff at ~ 3200 A . 
This will be useful in order to be consistent with what was 
done in the simulations. 

The one-dimensional distribution of Tea requires that 
the data be divided into broad bins (we will use five bins). 
In order to calculate Tea for each spectrum, we need to fit 
the continuum over the full Lya to Ly/3 region. Historically, 
the averaged value of the mean transmitted flux from many 
spectra has been used to calculate a global value of Tea (e.g., 
Kirkman et al. 2005). Fluctuations in Tea for entire spec¬ 
tra have not been estimated, however, presumably owing to 
problems with calculating a mean value on large scales where 
continuum fluctuations are seen to be comparable to fluctua¬ 
tions in the absorption (see for example the power spectrum 
of the raw spectrum in Hui et al. [2001] or McDonald et al. 
[2004]). The largest scale on which fluctuations intrinsic to 
the Lya forest have been measured is 153 Mpc at redshift 
2 = 1.9 by Tytler et al. (2004), who find that Lyman limit 
systems and metal lines make large contributions to fluctu¬ 
ations in the mean flux at this relatively low redshift. In our 
case, we will be working at 2 = 3, where these effects should 
be less important. Nevertheless, we expect that there will be 
a large contribution to the variation in Tea from sightline to 
sightline owing to continuum estimate errors as well as noise 
in the spectra and non-Lya absorption. We leave the ques¬ 
tion of how best to fit the continuum on the largest scales 
to other work (e.g., Suzuki et al. 2005), and in the mean¬ 
time, we will add a Gaussian distribution of errors to our 
simulation measurement to account for these effects when 
comparing our simulations with observations. The relatively 
large bin size should serve to mitigate effects of Tea errors 
in this first illustrative comparison. 

We fit the continuum using the often employed poly¬ 
nomial fitting method (e.g., Croft et al. 1998). A low order 
polynomial (3rd order in our case) is fit to the Lya forest 
region. We avoid the Lya and Ly/3 emission lines by re¬ 
stricting our fit to the rest wavelength range 1030-1190A. 
We then remove all points more than 2 a below the contin¬ 
uum and refit, iterating this procedure until convergence is 
achieved. Averaging over all Lya pixels in our dataset, we 
find a mean redshift 2 = 2.91. At this stage we divide out 
the continuum and compute a value of F. For all spectra, 
we find F = 0.8. Measurements of the mean flux and global 
value of Tea from higher resolution datasets better suited to 
this measurement find a lower value of F (e.g., Schaye et al. 
2004, Kirkman et al. 2005). Our measurement of F is on 
average biased high because we miss absorption in flatter 
low density regions around the continuum owing to our low 
resolution. We correct for this bias in a crude way by multi¬ 
plying the value of the fit continuum in each case by 1/0.85. 
This simple correction is adequate for our current purposes. 
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Figure 10. SDSS vs WMAP T ior z = 3. We show jackknife 
errors (10 subsamples). We have subtracted off a ring for back¬ 
ground from 15-25 arcmin around the pixel that the spectrum 
is in. The lines show predictions from LCDM after broadening 
the Tgff values with observational errors. The solid line is for the 
LCDM baryon density, dotted for LCDM X 5 (i.e. D 5 j = 0.4) 
and dash-dotted LCDM X 20 (i.e. Qbj = 0.8). 

having been set so that we reproduce the global Tefi value of 
Schaye et al. (2004) at a = 3. We make the assumption that 
any random variations in this bias from spectrum to spec¬ 
trum will simply act as an additional source of noise which 
we can model in our simulation data. 

In the bottom panel of Figure ITHl we have plotted a his¬ 
togram of the Teff values derived from the SDSS spectra. 
The mode of the distribution is =0.35, the same as the 
simulation results at « = 3 plotted in Figure |7| In the ob¬ 
servational data, there is a tail extending down to Tefr ~ 0.1 
and up to Teff ~ 1 which owes to noise and the continuum 
fitting uncertainties described above. 

4.3 CMB data 

We make use of the WMAP first year data (Bennett et al. 
2003) in our analysis. For the illustrative application of Lya 
-CMB cross-correlation in this paper, we employ the fore¬ 
ground cleaned version of the data presented by Tegmark 
et al. (2004). The CMB temperature is given in HEALPix 
format (12 x 512^) pixels, at the resolution of the WMAP 
W-band data (13 arcsec FWHM). 

The Tegmark et al. map is chosen by us because it 
combines the 5 WMAP channels into a map with the high¬ 
est angular resolution possible. However, because it is de¬ 


rived from channels with different frequencies and the ther¬ 
mal SZ effect has a weak frequency dependence, the map 
will have a scaled value of the SZ effect present, depen¬ 
dent on the weights used to combine the channels. In or¬ 
der to understand what impact this will have, we can use 
Equation 3 to examine the thermal SZ temperature fluc¬ 
tuations in the different bands. Equation 3 can be written 
in the form AT/T = f{i') x y, and at the lowest frequen¬ 
cies, f{i') —> —2.0. The five WMAP channels are centered 
at frequencies of 23, 33, 41, 61 and 94 GHz. The values of 
f{u) for the channels are then —1.97,—1.94,—1.91,—1.81 
and —1.56, respectively. The weights for the Tegmark map 
vary spatially and are not publically available, but one can 
see the maximum difference of the effect compared to our 
Rayleigh Jeans approximation {fO = —2.0) by considering 
the 94 GHz W-band map on its own, which gives approx¬ 
imately a 22% reduction in the predicted tSZ signal. If we 
instead assume that the weights will have a similar effect 
to those used by the WMAP team itself for their linearly 
combined map, we find by summing over the 5 bands that 
/effective = Wifi(u) = -1.79, or a ~ 10% reduction. The 
effect of frequency dependence is therefore small enough and 
uncertain enough that we neglect it in our illustrative ap¬ 
plication to observational data. In the future, however, as 
we have mentioned before, the frequency dependence can 
be used to better extract the signal. We will return to this 
point. 

For each pixel centered on one of the SDSS quasar spec¬ 
tra, we extract the corresponding CMB temperature of that 
pixel. In order to crudely mitigate the effects of primary 
CMB anisotropies, which act as noise for our signal, we sub¬ 
tract a local background consisting of the CMB temperature 
averaged over all pixels in an annulus of inner radius 15 ar¬ 
cmin and outer radius 25 arcmin, centered on the pixel in 
question. Our motivation for this background subtraction is 
the fact that the SZ signals fall off rapidly on large scales 
(e.g. see Figure]^, unlike the CMB primary signal. Because 
we will only be computing an upper limit to the SZ signal, 
this subtraction is a conservative assumption. 

The presence of CMB primordial anisotropies will pro¬ 
duce fluctuations on large scales. For example, there will be 
coherent fluctuations on the 100 arcmin scales which our 
simulated sky maps cover. The annulus subtraction tech¬ 
nique represents a crude way to deal with this. As we shall 
see below, the noise level on the pixel scale is sufficiently 
large with the WMAP data that we are not limited by the 
primary fluctuations in our analysis. However, in the future, 
with better data it will become imperative to model and re¬ 
move the effect of the primordial anisotropies. This could be 
done with matched filtering (e.g., Schaefer et al. 2004), and 
the technique tested with simulations of the CMB sky which 
include large-scale anisotropy, foregrounds, noise and an in¬ 
strumental beam response. We leave this for future work, 
and in Section El we will discuss other more sophisticated 
ways to deal with the CMB data. 

4.4 Results 

We find the rms temperature fluctuations in all 3488 pixels 
containing spectra to be 60/iK. In the middle panel of Figure 
HOI we show the rms temperature in each of 5 Tefi bins. The 
error bars range from 1.3 to 6.1/i K. This is close to the 













Lyman-alpha forest-CMB cross-correlation 11 


Poisson error estimate in each case, which is not surprising 
because the pixels are sampled quite sparsely from the data. 
Given that we expect ~ l/rK rms fluctuations at most from 
the kinetic SZ effect, it is clear that this particular dataset 
will fall short of giving us any meaningful constraint from the 
kSZ. Because the kSZ signal will add in quadrature to the 
rms from the noise and CMB primary signal, the highest Ten 
bin might have a total rms of \/60^ + l^/rK and the lowest 
\/602 — a difference of 0.02 /rK. In practice, because 

the CMB data is smoothed, we expect at least a factor of 10 
less signal (see the middle panel Figure |7J. The kSZ signal 
is therefore at least 500 times too small to be detectable in 
this way with these data. 

If instead we turn to the thermal SZ signal, the situa¬ 
tion is much more promising. As a mean temperature rather 
than a variance, it is easier to constrain this signal from noisy 
data. In the top panel of Figure ITT11 we show the mean tem¬ 
perature in pixels for the 5 bins of Lya Tefi . Owing to the 
background subtraction, the mean value of the pixels is zero. 
When comparing to the simulation data, we will renormalise 
so that the same is true. The observational error bars vary 
from 1.3 to 8.0 fiK per bin. 

From Figure]?! we expect the mean temperature decre¬ 
ment from the thermal SZ effect at redshift 2 = 3 to exhibit 
a total variation of ~ 0.6pK amongst the pixels with high¬ 
est and lowest Lya forest Tefi . This is for CMB maps at the 
WMAP resolution, but unfortunately not including noise in 
the Lya Tea estimates. As detailed in Section 4.2 above, 
we model the effects of noise in the spectra and continuum 
fitting uncertainties by adding a Gaussian-distributed error 
to the Tea from simulations. We add a Tea value randomly 
drawn from a normal distribution with standard deviation 
a = 0.1. Doing this to the simulated values of Tea enables 
us to reproduce the width of the histogram of observed Tefi 
values. In the top panel of Figure EOl we plot the curve of 
simulation CMB temperature (T) vs. simulated Tea after 
this noise has been added. The line which lies closest to the 
(T) = 0 axis is this prediction. As we can see, the predicted 
variation of (T) with Tea is barely noticeable on this scale, 
indicating that tight constraints lie only in the future with 
larger observational datasets. 

Nevertheless, the comparison is interesting because it 
constitutes an attempt to detect the ionised baryons at 
2 = 3, and the sensitivity increase required for success if 
our current cosmological models are correct is only about 
one order of magnitude. In order to place an upper limit on 
the temperature and density of ionised baryons using exist¬ 
ing data, one would ideally use a grid of simulations with 
different values of the baryon density, with the amplitude 
of matter fluctuations normalised so that their clustering 
and other properties are consistent with observations (e.g. 
Spergel et al. 2003, Tegmark et al. 2004). In this manner, 
both the clustering of baryons and the variation of mean 
IGM temperature with Dj, could be marginalised over to 
yield a constraint on Lib from ionised baryons, f2b,i. The 
effects of varying Qb on structure at high redshift were ex¬ 
plored with a suite of simulations by Gardner et al. (2003). 
These authors found for example that Lya forest absorption 
is virtually unaffected by changing Qb (by a factor ~ 6), as 
long as the ionising background intensity is adjusted to re¬ 
produce the same value of Tea ■ 

In the present case, because of the large error bars, the 


value of Db,i consistent with observations is likely to be so 
high that a simpler and more crude approach than a full 
simulation grid is all that is warranted. In order to model 
the effect of increasing Qb we simply multiply the simulation 
value of SZ (T) by a constant factor, assumed to be equiva¬ 
lent to multiplying Qb by the same factor. We do not adjust 
the Lya observables, imagining that any increase in baryon 
density can be cancelled out by the presence of a UV back¬ 
ground radiation field higher in intensity by the same factor. 
Moreover, we assume that increasing Qb does not affect the 
amplitude of matter fluctuations and that the temperature 
of the IGM stays the same as in the Qb = 0.04 simulation. 
A model with a very large value of Qb is likely to have many 
other differences related to gas physics and cooling, such as 
the number of collapsed baryonic objects. Additionally, of 
course, it will be incompatible with Big Bang Nucleosynthe¬ 
sis constraints (e.g.. Buries and Tytler 1998), and measure¬ 
ments of Qb from the CMB acoustic peaks (Spergel et al. 
2003). None of these problems concern us here, as our com¬ 
parison to simulations is essentially illustrative. 

Under these assumptions, a fit of the simulation data 
to the WMAP/SDSS data, both binned in the same manner, 
yields fib,: < 0.55 at 90% confidence and Db,i < 0.80 at 95% 
confidence. The limits on the fraction of the critical density 
in ionised baryons at 2 = 3 are therefore very weak. This 
is as one would expect from looking at Figure ESI where 
we show the simulation curves for various values of S2b,i- It 
is possible to conclude that the Universe is not closed by 
ionised baryons, but not much else. This situation is likely 
to be much improved with future observational data, as we 
will discuss below. 

Much of the statistical power in the measurement comes 
from the bin with the lowest value of Tea ■ Most of the pixels 
have moderate values of Tea but these pixels in the lowest bin 
also have low SZ decrement amplitudes. In this manner, the 
CMB-Lya forest cross-correlation can be said to be probing 
the low density IGM. 

5 SUMMARY AND DISCUSSION 
5.1 Summary 

Using a cosmological hydrodynamic simulation we have in¬ 
vestigated the CMB temperature fluctuations caused by the 
diffuse IGM at redshifts 2 = 2 — 6, also responsible for the 
Lya forest. Our main findings can be summarised as follows: 

(1) At redshift 2 = 2, the kinetic Sunyaev-Zel’dovich 
effect from gas within the Lya to Ly/3 region of the forest 
is predicted to cause CMB fluctuations of amplitude 
0.8 /rK rms averaged in a top-hat filter of angular radius 
1 arcmin. This remains roughly constant for higher redshifts. 

(2) The thermal SZ at redshift 2 = 2 yields an rms T 
of 2.4 /rK. At higher redshift, the relative fluctuation 
owing to the thermal SZ compared to kinetic SZ declines 
as expected because of the decreasing mean temperature. 
However, there is still a substantial Trms at higher redshift, 
0.4 nK at 2 = 4. 

(3) There is a correlation between Lya forest proper¬ 
ties of quasar spectra with CMB temperature fluctuations 




12 R.A.C. Croft et al. 


caused by both the kinetic and thermal SZ effects. In 
particular, there is a strong relation between the mean 
effective optical depth reft averaged over the entire Lya 
forest region and the CMB temperature. Pixels (of size ~ 1 
arcmin) in the top 10 percentile of reft are predicted to have 
a mean CMB temperature 1.5 fiK lower than the bottom 
10 percentile at a = 2, owing to the thermal SZ effect. The 
quantity reft is relatively sensitive to density fluctuations in 
the diffuse IGM. 

(4) There is, however, little correlation between the 
total HI column density in spectra and the CMB temper¬ 
ature, with no relation at all seen for column densities 
Nhi > 10^®cm“^. This statistic is most sensitive to the 
neutral hydrogen density in saturated regions, indicating 
that SZ effects correlate best with relatively diffuse unsatu¬ 
rated gas. 

(5) We have compared the correlation of Lya Tefr with 
CMB T seen in simulations of the ACDM cosmology with 
observational data from the Sloan Digital Sky Survey and 
WMAP. Because the kinetic SZ signal and pixel noise add 
in quadrature, no constraints are possible on the IGM 
velocity field with the present data. However, the sensitivity 
is within an order of magnitude of that necessary to detect 
the predicted thermal SZ effect from the Lya forest. We 
are able to place a weak limit on the density in ionised 
baryons at 2 = 3 of Hb,! < 0.80 at 95% confidence. Below 
we outline future developments which would make a direct 
detection possible. 

5.2 Discussion 

In the present paper, we have examined the relationship 
between Lya forest absorption and the SZ effect in general, 
without the expectation that a detection of a correlation can 
be made from current data. In the future, the quality and 
size of datasets will increase, and it will become worthwhile 
to plan a strategy for the cross-correlation that makes the 
best use of the data. Two obvious improvements could be 
made in the future. In our illustrative comparison of WMAP 
and SDSS data we did not use the full angular or frequency 
information, both of which would help with the sensitivity 
of the test. 

As we have seen that the thermal SZ effect offers the 
best hope for detection of the ionised baryons, making use 
of the spectral difference between the SZ decrement and the 
primary CMB signal will enable easier detection of the ef¬ 
fect. One could imagine fitting the Tefr -T correlation in dif¬ 
ferent CMB frequency bands to do this. For example, in the 
work of Huffenberger et al. (2004), a method was derived 
for combining the WMAP channels in an optimal manner 
to remove the primary anisotropy, effects of noise and con¬ 
taminating point sources, enabling the thermal SZ signature 
to be constrained directly. In this paper we have concerned 
ourselves with predictions for the Rayleigh Jeans part of the 
spectrum, which is reasonable given the large observational 
uncertainties in present data. We note that at higher fre¬ 
quencies the thermal SZ signal becomes positive (for exam¬ 
ple at 330 GHz the magnitude is the same as the Rayleigh 
Jeans signal but with the opposite sign). Bolometers tend 
to be more sensitive than radiometers at present, and they 


operate at frequencies > 100 GHz. These differences could 
be used profitably to increase the signal to noise of any de¬ 
tection. Also, with the WMAP data, we have limited our 
analysis to single pixels. Making use of the full spatial cross¬ 
correlation in our detection of the signal will be useful when 
CMB datasets with higher angular resolution become avail¬ 
able. Another possibility could be to use wavelets to isolate 
the relevant physical scales in a future analysis (see e.g., 
Vielva et al. 2005). 

One can ask whether emission from the QSOs them¬ 
selves might create a signal, as the SED of QSOs in the 
microwave band is not well constrained (Perna & Di Mat- 
teo 2000, White & Majumdar 2004). Apart from the pos¬ 
sible intrinsic microwave emission of QSOs, there will also 
be an intrinsic thermal SZ effect from galaxy clusters which 
host QSOs. Neither of these will necessarily bias any mea¬ 
surement of the Lya -CMB correlation unless QSOs with 
more Lya absorption in their spectra have more intrinsic 
microwave emission. This seems unlikely in the former case, 
although in the case of an SZ signal associated with clus¬ 
ters in which the QSOs lie one could argue that this is part 
of the signal that we are trying to measure. Otherwise, it 
could be removed if we so desire by not using the part of the 
spectrum closest to the QSO. 

The method we have outlined in this paper is one way 
to potentially detect ionised baryons in the IGM at high 
redshift. In general, one would like to know if there are any 
other methods that can be used. At lower redshifts, for ex¬ 
ample (2 < 1), the ionised baryons are mainly collisionally 
ionised, with such a low neutral fraction that they cannot be 
seen in the Lya forest. Their higher temperature, however, 
means that this Warm Hot Intergalactic Medium (WHIM, 
see e.g., Cen et al. 1995, 1999, Dave et al 2001) can po¬ 
tentially be detected by looking for X-ray absorption (e.g.. 
Fang et al. 2002) or emission (Yoshikawa et al. 2004, Fang 
et al. 2005) or perhaps in Lya emission (e.g. Furlanetto et 
al. 2003, 2005). The thermal SZ signal at low redshifts is 
also one way to detect the WHIM. At higher redshifts and 
lower IGM temperatures considered in this paper however, 
the X-ray signal (see e.g., Croft et al. 2001) is not likely to 
be observable. 

The overall y-distortion in the CMB spectrum from 
the hot IGM should eventually be detectable. Limits from 
the COBE FIRAS instrument (e.g., Mather et al. 1990) are 
y < 2.5 X 10“® at 95% confidence, more than an order of 
magnitude above the mean y caused by the high -2 IGM. 
However, separating out the contribution from the high 2 
IGM will be difficult without a spatial correlation of the 
type we have proposed here. One could also use other trac¬ 
ers of the high 2 density field for a cross-correlation, such as 
Lyman-break galaxies, or QSOs themselves, although their 
connection to the IGM would be less easy to interpret. 

Another interesting question to ask is how likely it is 
that we will be able to constrain the high 2 IGM velocity 
field, given that the kinetic SZ signal from the photoionised 
IGM is expected to be relatively important (Loeb 1996). 
Because of the possible -|-ve or -ve nature of the effect, de¬ 
pending on the velocity, we cannot use a mean signal, unlike 
the thermal SZ. By using an rms fluctuation as our statistic, 
we therefore need a noise level which is extremely low. Since 
the spectrum of the kSZ is the same as the CMB, the pri¬ 
mary fluctuations become even more important sources of 
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noise than with the thermal effect. Given that the rms fluc¬ 
tuations on scales of 1 arcmin from the kSZ are predicted 
to be ~ 5/iK at 2 = 2, one would need the “noise” from 
primary anisotropies and other sources to be less than this 
in order to detect a signal. 

The thermal SZ is more promising. Future experiments 
which will be useful for improving constraints on in¬ 
clude the Planck satellite mission^, the Sunyaev Zel’dovich 
Array (SZA®), the Atacama Cosmology Telescope (ACT, 
see Kosowsky 2003), and the South Pole Telescope (SPT'*). 
Given the expected noise characteristics and angular resolu¬ 
tion/sky coverage of these instruments, we expect substan¬ 
tial improvements over the constraints derived using WMAP 
in this paper. We have seen in the WMAP case (Section 
that we are limited both by the noise level (~ 60/rK rms) 
and the resolution of the data. As an example, the ACT 
angular resolution is projected to be ~ 1 arcmin with an 
rms detector noise level of ~ 2/rK per pixel in the 145 GHz 
band. Even without allowing for the higher resolution or us¬ 
ing the full angular cross-correlation information, one would 
therefore expect to extract the same level Lya forest ther¬ 
mal SZ signal using (60/2)^ fewer spectra (assuming the 
primary anisotropy can be filtered out in frequency space). 
The ACT signal will be lower by about 50% at 145 GHz 
than the Rayleigh Jeans prediction, but the other frequency 
bands can be used to increase the significance of the detec¬ 
tion. In order to yield a 2-a detection of Hb, using pixel-pdf 
statistics, as in this paper, we estimate that ~ 300 quasar 
Lyof forest spectra would be needed in the ACT survey area, 
assuming the BBN value of Qbh^ (0.019; Buries & Tytler 
1998) Using angular cross-correlation information would in¬ 
crease the significance of the detection. The lack of overlap 
between the ACT fields and the SDSS would necessitate a 
southern quasar survey, although it should be borne in mind 
that only extremely low resolution spectra are needed. 

We note that in our analysis we have assumed that the 
ratio of neutral to ionised hydrogen is that expected if the 
ionising background radiation field does not vary spatially 
in intensity. If there are coherent spatial fluctuations in the 
photoionisation rate (P in Equation 1) then the correlation 
between Lya optical depth and CMB temperature decre¬ 
ment will be affected by an additional source of scatter. 
This is likely to be small, however with the spatial cross¬ 
correlation reduced by 5% or less (see e.g., Croft 2004) and 
could in principle be calibrated out using theoretical mod¬ 
elling in the manner of McDonald et al. 2005. 

Finally, one can ask why it is important to try to detect 
the ionised component of the baryons at high redshift when 
the neutral hydrogen tracer has already been well measured 
(e.g., Rauch et al. 1997). One reason is that a detection 
would act as a consistency check and enable us to verify di¬ 
rectly the baryon inventory at these redshifts. Also, at the 
moment, the best constraints on the baryon density in ob¬ 
served structures (if we set aside for the moment relatively 
indirect CMB and Big Bang Nucleosynthesis, BBN measure¬ 
ments e.g., Spergel et al. 2003, Walker et al. 1991) comes 
from looking at the mean absorption level in the Lya for- 


^ http://planck.esa.int 
® http://astro.uchicago.edu/sza 
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est (e.g., Rauch et al. 1997, Weinberg et al. 1997, Hui et al. 
2002, Tytler et al. 2004). However, these measurements are 
degenerate with constraints on the ionising background ra¬ 
diation intensity J. Often, Hb from BBN or CMB is assumed 
and then the Lya forest measurement is used to constrain 
the ionising radiation intensity (e.g., Kirkman et al. 2005). 
If we use actual measurements of J (from e.g., the prox¬ 
imity effect e.g., Scott et al. 2000) then constraints on Hb 
directly seen in the Lya forest become weaker (J is only 
known to within roughly a factor of two). Measurements of 
the ionised baryon fraction such as those we propose may 
therefore become competitive. 
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